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ABSTRACT 



Aims. We use different tracers, such as Ha, ultraviolet (UV), and infrared (IR) emissions at various wavelengths, to study the dust and 
star-formation (SF) conditions throughout the disk of M 33. 

Methods. We derive the radial distribution of dust, of the old and young stellar population using Spitzer and GALEX data, comple- 
mented by ground-based optical data and available surveys of atomic and molecular gas. We separate the contribution of discrete 
sources to the IR brightness from the diffuse emission. 

Results. At 8 and 24 pm, discrete sources account for > 40% of the IR emission in the innermost 3 kpc, and for < 20% further 
out. We find that stochastic emission from very small grains in the diffuse interstellar medium accounts for only ~10% of the diffuse 
24 fim emission, and that dusty circumstellar shells of unresolved, evolved AGB stars (carbon stars) are a viable alternative. The 8 pm 
profile suggests that PAH emission declines faster with radius than the dust continuum. In annular regions, 0.24 kpc wide, we find a 
mean extinction value for stellar continuum Ay ~ 0.25 mag with a weak dependence on radius, consistent with the shallow metallicity 
gradient observed. Dust opacity derived from the 160 /jm emission decreases instead by a factor 10 from the center to edge of the SF 
disk. 

Conclusions. Using extinction corrected UV and Ha maps we find the global SF rate in M 33, over the last 100 Myr, to be 0.45 ±0.10 
M Q yr -1 . FIR (far-IR) and TIR (total-IR) luminosities can trace SF even though a high conversion factor is required to recover the 
effective rate. If carbon stars are powering the diffuse 24 pm emission in M 33 this can trace star formation 1 Gyr ago and provide a 
more complete view of the SF history of the galaxy. Today the SF rate declines radially with a scale length of ~ 2 kpc, longer than for 
the old stellar population, suggesting an inside-out growth of the disk. 

Key words. Galaxies: individual (M33) - Galaxies: ISM - Galaxies: Local Group - Galaxies: spiral 



1. Introduction 

Our knowledge of the conditions in the interstellar medium 
(ISM) that favour star formation (SF) is mostly based on our 
Galaxy and on luminous, gas-rich galaxies with high SF rates 
(SFRs). Thus, any proposed SF scenario is biased against dif- 
ferent Hubble types, especially against quiescent galaxies with a 
low molecular-gas content. Moreover, our non panoramic view 
of the Milky Way, with our solar system deeply embedded in it, 
does not allow us to determine firmly important properties of our 
host galaxy, such as the global SFR and its variations across the 
galactic disk. High sensitivity and resolution in external galax- 
ies are desirable in order to link properties integrated over the 
disk to the distribution of individual star-forming regions, to the 
ISM components and to large scale perturbations in the disk. 
Data from recent space missions, such as the Galaxy Evolution 
Explorer ( GALEX) and Spitzer satellite, allow us today to inves- 
tigate these properties in the nearest galaxies, before the next 
generation of telescopes will resolve star-forming sites in more 
distant ones. 

The Local Group galaxy M 33 has a high SFR per unit area 
compared to M 3 1 and a low extin ction towards SF regions, ow- 
ing to its mod erate gas content ( Corbelli 2003) and its low in- 
clination (54°, McConnac hie et alj2 006). It is a rather quiescent 
blue disk galaxy, showing no sign of rec ent mergers or interac- 
tion. At a distance D of only 840 kpc (see lFreedman et aT]|l99ll) 
it allows us to examine its disk with detailed resolution and con- 



stitutes a local prototype for studying the ISM and its relation to 
SF in blue, low luminosity objects. Observations of the molec- 
ular, atomic and ionis ed gas in M 33 have been carried out in 
the past few years (e.g ICorbelli & Saluccill200fA ICorbelli|[2003l: 
lEngargiola et alj200 3; Heyer et al. 2004) to gether with metallic - 
ity surveys (Magri ni et al.ll2007allbtlRosolowskv & S imon 2008) 
of the young and old stellar populations across its disk.The 
galaxy has about 10% of the total gas content in molecular form 
and a very shallow metallicity gradient from 0.5 kpc to the edge 
of the star-forming disk. 

In the first pap er of this series, dedicated to the SF in M 33 
dVerlev et al.ll2.Q07l hereafter Paper I), we used Spitzer maps to 
select a large sample of infrared (IR) emitting sources in the 
24 pm map of this galaxy and to identify their counterparts in 
all IRAC wavebands and in Ha. There is large spread in the cor- 
relation between the Ha and the IR fluxes of star-forming sites, 
similar to what has been found by Spitzer in the fai nt, metal de- 
ficient local dwarf NGC6822 dCannon et al.ll2006l) . In selected 
areas of this galaxy the Ha-to-IR ratios show variations as high 
as a factor of 10, co ntrary to the tight correlation found in th e 
bright galaxy, M 51 dCalzetti et alj|2005b iKennicutt et al.ll2007l) . 
In a forthcoming paper (Corbelli et al. 2008, in preparation) we 
will use multiwavelength photometry to characterise how stars 
form in individual star-forming sites in M 33. 

In this paper, the second of this series, we focus on the in- 
tegrated IR emission and on its radial variations to establish the 
role of discrete star-forming regions in powering the dust emis- 
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sion across the disk. Indeed, there has been a long debate on the 
role of localised warm sources to power the dust emission and 
the IR fl uxes in general. C om paring IRAS and Ha o bservations 
of M 33. lRice et afldl990l) and lDevereux et al.l(11997l) found that 
more than 70% of the far infrared (FIR) is linked to massive 
stars complexes . On th e other hand, based on ISO observations, 
Hipp elein et al.l (120031) suggested that a substantial part of the 
FIR emission could be due to the cold component heated by the 
diffuse interstellar radiation field (ISRF). Part of the difficulties 
of previous surveys in determining the grain heating source had 
to do with the difficulties in resolving sources from IRAS and 
ISO observations. With Spitzer, we can now probe very faint H n 
regions and tackle this problem more effectively. In Paper I we 
have shown that discrete sources account for less than half of the 
IR emission at 8 and 24 /mi. Understanding the origin of these 
high diffuse fractions in M 33 will be one of the goals of this 
paper. We will examine dust heating sources and the contribu- 
tion of dust in different environments (diffuse ISM, H n regions, 
evolved stars) to the IR emission at several wavelengths as we 
move from the center to the outermost regions of M 33. We then 
use M 33 as a prototype to test the reliability of UV continuum 
fluxes, recomb ination lines and IR continuum fluxes (see, e.g., 
iKennicutll 19981 for a complete review), as tracers of SF to pos- 
sibly understand more distant galaxies, similar to M 33 but un- 
resolved by present instruments. 

The plan of the paper is the following: in Sect.|2]we present 
the data, in Sect.[3]the integrated luminosities in UV, Ha and IR 
and their radial profiles, as well as the atomic and molecular gas 
distributions. Optical and UV extinction and dust abundance are 
presented in Sect. |4] In Sect. [5] we discuss the origin of the IR 
emission in discrete sources and in the more diffuse component, 
with particular emphasis to the 24 pm band. In Sect. [6] we dis- 
cuss different tracers of SFRs and Sect. [7] summarises our main 
results. 



2.3. Infrared data 

The MIR and FIR data of M33 are those provided by the 
Spitzer InfraRed Array Camera (IRAC) and Multiband Imaging 
Photometer for Spitzer (MIPS) ins truments (W erner et al.l l2004: 
iFazio et alJfeOO^ lRieke et al.ll2004l) . The complete set of IRAC 
(3.6, 4.5, 5.8, and 8.0/im) and MIPS ( 24, 70, and 160 /an) 
images of M33 is describ ed in iPaper It the Mopex software 
( Mako voz & Marleaull2005h was used to gather and reduce the 
Basic Calibrated Data (BCD). The complete field of view ob- 
served by Spitzer is very large and allows us to achieve a com- 
plete view of the star-forming disk of M 33, in spite of its rela- 
tively large extension on t he sky. T he images were background 
subtracted, as explained in IPaper il . The spatial resolution mea- 
sured on the images are 2'.'5, 2'.'9, 3'.'0, 3'.'0, 6", 16", and 40" 
for IRAC 3.6, 4.5, 5.8, 8.0 Aim and MIPS 24, 70, and 160 fim, 
respectively. The absolute photom etric calibration uncertainties 
are be tter than 10% for all IRAC dFazio et al.ll2004tlReach et all 
120051) and the MIPS 24 pm bands (Engelbr acht et al.ll2007h and 
are better than 2 0% for the MIPS 70 and 160/im channels 
dRieke et alj|2004 . 

2.4. 21 -cm and millimeter data 

The radio and mm databases to infer the atomic and 
molecular gas distributions are the f ollowing: Westerbork 
Radio Synthesis Telescope (WRST) array dDeul & van der Hulstl 
1987, 24" x 48" spatial re solution) and Arecibo single dish 
dCorbelli & Schneidel[l997l 4' sp. res.) 21-cm line data for the 
atomi c gas; Berkeley Illino is Maryland Association (BIMA) ar- 
ray dEngargiola et al.ll2.Q03L 13" sp. res.) and Fiv e College Radio 
Astronomy Observatory (FCRAO) sing le dish dCorbellill200l 
45" sp. res.) for molecular gas emission as traced by the CO J=l- 
rotational line. 



2. Observations and data sets 



2.1. Ultraviolet data 



To investigate the UV continuum emission of M 33, we make 
use of ob servations by the G alaxy Evolution Explorer ( GALEX) 
mission dMartin et all I2005I). We use the data distributed by 
iGil de Paz et al. ( 2007b . A description of GALEX observations 
of M 33 in the far-UV (FUV, 135 0-1750 A) and near-U V (NUV, 
1750-2750 A) can be found in iThilker et all d2005l) . together 
with the relevant information on data reduction and calibration 
procedures. 



2.2. Ha data 



3. Integrated emission and radial profiles 

The s udden drop of the Ha emission observed by iKennicutj 
(1989) in M 33 , and shown also here later in this Section, de- 
fines the extent of the star-forming disk to be 7 kpc in radius. In 
this Section, we first estimate the total luminosity of M 33 in UV, 
Ha and IR wavelengths by integrating the observed fluxes over 
the star-forming disk; next we extract radial profiles of ellipti- 
cally averaged fluxes at IRAC and MIPS wavelengths and also 
from the FUV, NUV, and Ha images. We will use these profiles 
to evaluate the radial dependence of the average extinction which 
we then use to derive extinction corrected profiles and total lu- 
minosities in UV and Ha. The IR emission (IRAC+MIPS) will 
be considered optically thin: luminosities are not corrected, and 
face-on values for the surface brightness are given by correcting 
only for the disk inclination to the line of sight (54°). 



We use th e narro w-band Ha image of M33 obtained by 
iGreenawalj {[998). The reduction procedure, using standard 
IRAF3 tasks to subtra ct the con tinuum emission, is described 
in detail in iHoopes & Walterbosl (l2000h . The total field of view 
of the image is 1.75 x 1.75 deg 2 (2048 x 2048 pixels of 2'.'0). 



3.1. Integrated luminosities 

The IR luminosities within 7 kpc for the different IRAC and 
MIPS bandpassefl are given in Table [T] together with the 3- 
1000 //m total- i nfrar ed luminosity (TIR) derived according to the 
iDale & Heloul d2002) relation: 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



F(TIR) = 10~ 14 x (19.5 F v (24) + 3.3 F v (70) + 2.6 F„(160)), (1) 



- Following ICalzetti et all J2005t) , the convention 
F, I, L(band) = vF v , vl v , vL„(band) is adopted. 
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Table 1. Total fluxes and luminosities within 7 kpc. Ha, FUV 
and NUV luminosities corrected for extinction are also shown. 



Wavelength 


Total fluxes 


T 

L, 


L(corr) 


[fim] 


LJyJ 


[erg s" 1 ] 


[erg s~'] 


3.6 


19.9 


1.4 10 




4.3 


1 o o 

iz.o 


1.1 iU 




J.O 


24.4 


1.1 10 42 




8.0 


51.8 


1.6 10 42 




24 


49.4 


5.2 10 41 




70 


666. 


2.4 10 42 




160 


1909. 


3.0 10 42 




TIR 




6.8 10 42 




Ho- 




3.1 10 40 


4.2 10 40 


FUV 




3.9 10 42 


7.0 10 42 


NUV 




3.2 10 42 


6.3 10 42 



with F(TIR) in W m 2 and F v in Jy. The TIR luminosity of 
M33 is 1.8 x 1 9 L . T he TIR luminosity we estimated from 
8 and 24 pm in iPaper 5 was 1.0 x 10 9 L Q , but within a shorter 
radius (5 kpc), hence the two values are quite consistent. Using 
the values in Table Q] we show in Fig. Q] the mid- and far-IR 
spectral energy d istribution (SED) together with IRAS data from 
iRice et alJ (I 19901) which, given the uncertainties, turn out to be 
in satisfactory agreement with our measurements. In the lower 
panel of Fig.Q] the SED is depicted as vl v which can be directly 
comp ared with the IR SED of 71 nearby galaxies by lDale et al.l 
(120051) . The SED results from two major components: the stel- 
lar photospheric emission, which decreases from 3.6 to longer 
wavelengths, and the rising continuum due to dust emission. The 
peak at 8.0/itn gives clear evidence of the presence of PAHs in 
M33. 

The integrated FUV and NUV luminosities within 7 kpc 
are 3.9 x 10 42 e rg s" 1 and 3.2 x 10 42 erg s _1 respectively 
dGil de Paz et al.l 120071). T he Ha luminosity (as given by 
iHoopes & Walterbosll2000l) is 3.1 x 10 4() erg s 1 . Again, TableQ] 
summarises the total fluxes, and luminosities measured within 
the SF disk of M33 before applying any extinction correction, 
and with extinction corrections derived in Sect. 14. II 

3.2. IR radial profiles 

Previous work has suggested that the SF disk in M33 shows 
three morphologically distinct areas: first, the innermost 1 kpc 
where a weak bar is in place and metallicity is enhanced; second, 
the inner 3.5 kpc disk where spiral arms and GMCs are promi- 
nent; and finally the outer disk where structures weaken before 
the edge of the SF disk occurs at about 7 kpc. It is therefore of in- 
terest to examine how the old and young stellar populations and 
the dust emission vary among these three regions and in general 
as a function of radius. 

To extract the IR radial profiles from images, we used the 
IRAF task ellipse, fixing the center to the nominal position of 
M 33 (01 h 33 m 50 s .90, +30°39'35'/8), the axial ratio to 0.59 as for 
a fixed disk inclination of 54° and the position angle to 22.5°. 

We show in Fig. [2] the elliptically averaged surface bright- 
ness (F0 in the IRAC and MIPS bands, as a function of galac- 
tocentric radius after correction flux densities to face on values. 
In Fig. |2] radial bins increase logarithmically: the multiplicative 
factor between two consecutive elliptical semiaxes is 1.1, that is 
an increment of 10%. The errors are computed by quadratically 
adding the error of the mean along the ellipse provided by el- 
lipse task itself and the error on the subtracted background. This 




Wavelength [pm] 




Wavelength [/jm] 



Fig. 1. The infrared spectral energy distribution of M 33 is 
depicted by (black) pluses (IRAC and MIPS Spitzer data). 
The IRAS data (12, 25, 60, and 100 ^m, green circles) from 
IRice et ail (1 19901) are compatible with our data within the quoted 
uncertainties. 



background error is estimated for each image by considering the 
noise and flatness of the image in the most peripheral areas (con- 
sidering 30 different areas around the borders). These estimates 
of the background uncertainty agree pretty well with the level of 
oscillation of the mean intensity from ellipse when it achieves 
background flux levels. 

The extraction of the IRAC and MIPS-24 profiles was per- 
formed on images where foreground stars had been previously 
removed; we observed, in fact, that they introduce significant lo- 
cal enhancement in the profiles, particularly at 3.6 and 4.5pm, 
the most sensitive to photospheric emission. At all wavelengths, 
the average surface brightness shows an enhancement in the 
innermost 0.5 kpc region (see Sect. [33}, then declines out- 
wards exponentially. While the old stellar population dominates 
at short wavelengths, PAHs and Very Small Grains (VSGs) in 
the surrounding of hot stars are expected to contribute at 8 and 
24 pm respectively. Hence these last profiles display more fea- 
tures, especially in the inner disk where spiral arms are promi- 
nent. These are visible also in the emission at 70 pm, while at 
160/mi the coarse resolution limits the detection of such fea- 
tures in the elliptical averages. 
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Fig. 2. Radial profiles, using logarithmical bins, of the Spitzer 
MIPS+IRAC surface brightness in arbitrary units. For clarity, re- 
sults for the different bandpasses have been shifted in the vertical 
direction (within brackets, the logarithmic shift to be added to 
the ordinates to obtain the real surface brightness in mJy pc~ 2 ). 
The dashed vertical line at 7 kpc marks the edge of the star- 
forming disk, see text. 

Table 2. Scale lengths in kpc for the radially declining surface 
brightness at various IR wavelengths. 



Wavelength 


0.5 : 


3.5 kpc 


3.5 : 


7 


kpc 


0.5 : 


7 kpc 


3.6 pm 


1.51 


±0.02 


1.61 


± 


0.05 


1.56 


±0.01 


4.5 fim 


1.52 


±0.04 


1.60 


± 


0.15 


1.55 


±0.03 


5.8 yum 


1.59 


±0.06 


1.36 


± 


0.13 


1.61 


±0.04 


8.0yum 


1.52 


±0.07 


1.21 


± 


0.05 


1.44 


±0.02 


24 pm 


1.71 


±0.10 


1.59 


± 


0.20 


1.77 


±0.06 


70 pm 


1.78 


±0.10 


1.55 


± 


0.14 


1.74 


±0.05 


160 pm 


2.41 


±0.09 


1.48 


± 


0.04 


1.99 


±0.02 



The x 1 fit exponential scale lengths fitted from 0.5 kpc to 
the edge of the star-forming disk are reported in Table [2] The 
scale lengths are computed after rebbining the data into ra- 
dial bins of equal size (50 pixels i.e. 0.24 kpc) taking into 
account the relative brightness uncertainties). All MIPS scale 
lengths are, within 3 cr compatible with those determined by 
Tabatabae i et al.l d2007b) (see their Table 7). The profiles at 3.6, 
4.5, 24, and 70 //m are consistent with a single exponential scale 
length throughout the star-forming disk. This is ~ 1 .55 kpc for the 
old stellar population emission and ~ 1 .75 kpc for the dust emis- 
sion at 24 and 70 pm. Two different slopes for the inner (0.5- 
3.5 kpc) and outer (3.5-7.0 kpc) star-forming disk are found in- 
stead at 5.8, 8.0 and 160pm. At these wavelengths the IR bright- 
ness declines more steeply beyond 3.5 kpc and drops at the edge 
of the star-forming disk. Drops at the disk edge are less evident 
at other wavelengths, even though the uncertainties at large radii 
make the data consistent either with the no-drop hypothesis or 
with the one that undoubtedly exists, namely at 160 pm. We also 
note that the disk of M 33 starts warping beyond 7 kpc, which 
adds uncertainties which have not been taken into account. 

The similarity of the 5.8, 8.0, and 160pm emission may be 
understood by comparing radial trends of PAH emission and 
cold dust. PAHs contribute to the 5.8 un\ emissi on a nd even more 
tothe 8.0pm emission band. lHaas et al.l d2002l) and lBendo et al.l 
(2008) have shown that PAH and cold dust emission correlate 
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Fig. 3. Brightness radial profiles in the FUV (small crosses), 
NUV (small pluses) and Ha (small diamonds). Extinction cor- 
rected brightness profiles in bins of 0.24 kpc are also shown in 
the FUV (large crosses), NUV (large pluses) and Ho- (large di- 
amonds). Ft is here converted to an absolute surface luminosity 
scale, that is multiplied by v x AnD 2 . The 8 pm radial profile 
(squares) is plotted for comparison with the UV brightness. 



spatially and inferred that the diffuse ISRF is the power source 
of both emissions. Indeed, M33 both 8pm and 160pm show 
a steeper decline in the outer disk. This might be due to a de- 
crease in the surface density of the 8 and 160 pm carriers since 
the decrease of the energy density of the ISRF in the outer disk 
is not steeper (see also Sec. I3.31 >. However, the 8pm profile de- 
clines with a steeper slope throughout the whole disk implying 
that the 8 to 160pm ratio also declines as we move outwards. 
One possibility is that the PAH abundance or emission is more 
sensitive to quantities which vary across the disk, such as the 
metallicity, the efficiency of supern ovae shocks in destroying the 
molecules dO'Halloran et al.ll2006l) . or the disk gravity, which 
affects the cooling and ionisation balance in the ISM (see Sect. 
13.4b . M33 has a shallow metallicity gradient, and the oxygen 
abundance [O/H] is always above -3.8, except at the edge of the 
star-forming disk (Mag rrnT et al.ll200"7bT) . This casts some doubt 
that the metallicity is the main driver of the fast radial decline of 
PAH emission. 

The scale length inferred for the old stellar population us- 
ing the 3.6 and 4.5pm emission is consistent with what earlier 
studies have derived in the inner part of the disk usin g photom- 
etry in the K-band (1.4+0.1 kpc. lRegan & Vogefl994l) . Ground 
based near IR data for this galaxy does not extend radially out- 
wards as far as Spitzer data, but a colour ma p of the 3. 6-2. 2pm 
emission using the 2MASS (J+H+K) image dJarrett et alj|2003l) 
reveals a radial gradient between 0.5 and 5 kpc due to the smaller 
scale length of the J+H+K emission (of order 1 kpc). This might 
imply that intermediate age, cool supergiants contribute a sub- 
stantial fraction of the NIR emission. The colour map follows 
the spiral morphology, suggesting that cool supergiants are more 
likely candidates than giants. 

3.3. UV and Ha line radial profiles 

To evaluate the recent SFR across M 33 we will use the FUV and 
Ha radial profiles. It is more reliable to compute the SFR from 
integrated emission across the disk rather than from single SF 
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Table 3. Scale lengths for the Ha, UV (uncorrected and cor- 
rected for extinction), the molecular and total gas emission. 



Line 


Wavelength 


Range 


K 


/i,(corr) 






kpc 


kpc 


kpc 


Ha 


6563 A 


[0.5 : 6.0] 


2.0+0.2 


1.8+0.1 


NUV 


1750-2750 A 


[ 0.5 : 6.0 ] 


2.3+0.1 


2.1+0.1 


FUV 


1350-1750 A 


[ 0.5 : 6.0 ] 


2.4+0.1 


2.2+0.1 


CO BIMA 




[0.5 : 6.0] 


1.4+0.1 




CO FCRAO 




[ 0.5 : 6.0 ] 


2 4 +1 




Total gas 




[0.5 : 6.0] 


16!> 8 





sites. This is because SF sites lose part of their energetic photons 
in the surrounding ISM and the stars themselves diffuse out of 
their birthplace over time. FUV is a diagnostic of a somewhat 
younger (30-100 Myr) star formation than NUV (100-300 Myr), 
while Ha relates to the youngest stellar population. In Fig. [3] 
we show t he observed radial profiles of the NUV, FUV and Ha 
emission dThilker et al.ll2005l) in bins of 0.24 kpc. Ha, FUV and 
NUV display very similar radial profiles which points to a scarce 
radial ev olution of the SFR i n the last 300 Myr. Unlike other 
galaxies (Thil ker et alj|2007ab . in M33, the UV emission does 
not extend farther than Ha. In addition, we observe no strong 
gradient in the colour ratio FUV-NUV; M 33 is only slightly 
bluer with incre asing radius opposit e to the more pronounced 
trends found by iBianchi etall (120051) in M 51 and M 101. This 
confirms the lack of an appreciable age gradient in the young 
stellar populations, and precludes strong radial variations of the 
UV spectrum. 

UV and Ha show a strong enhancement in the innermost 
0.5 kpc likely d ue to the effect of the bar fu eling of the circum- 
nuclear region (ICorbelli & Walterbosll2007l) and that is another 
reason why we excluded that region from our global fits. The ra- 
dial scale lengths of the UV and Ha brightness are shown in the 
fourth column of Table 3. They do not have significant variations 
between the inner and the outer disk and hence we display only 
their average values across the whole disk. In the last column 
of Table 3 we give these value after correcting for extinction, as 
described in the next Section. 

All three young SF tracers -Ha, FUV, and NUV- display 
a compatible scale length of about 2 kpc while the oldest stel- 
lar populations (3.6 and 4.5/mi) exhibit a shorter one (approx- 
imately 1.5 kpc). The fact that the youngest stellar component 
has the longest scale len gth supports the insid e-out disk growth 
scenario for M33 (e.g. Mag rini et al.l l2007ab . In general, our 
data suggest that dust emission at MIPS wavelengths in M 33 
is also declining radially with longer scale lengths than the old 
stellar population, especially in the inner disk. The average dust 
emission scale length over the SF disk is consistent with the one 
traced by the young stellar populations even though the shorter 
scale length of the IR in the outer disk, especially at 160 /mi 
suggests a decline of the dust surface density. 

3.4. Atomic and molecular gas distributions 

In Fig. |4] we show the atomic and molecular gas column den- 
sity as a function of galactocentric radius and Table 3 gives 
their radial scale lengths. The total gas distribution, sum of the 
H i (using the Arecibo data, except for the first four inner radii 
where we used the values from Westerbork) and of the diffuse 
molecular gas (using the FCRAO data), is also shown in Fig. [4] 
Uncertainties on the radial scale lengths for the CO-FCRAO 
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Fig. 4. Gas surface densities. The atomic gas data is shown by the 
open squares and open triangles for the Arecibo and Westerbork 
data, respectively. The molecular gas is shown by crosses and 
pluses for the BIMA and FCRAO data, respectively. The sum of 
the atomic (using the Arecibo data, except for the first four inner 
radii where we used the values from Westerbork) and molecu- 
lar gas (FCRAO) is also displayed with filled circles. Error bars 
show the dispersion in each bin. 

and total gas are computed using the x 1 statistic taking into ac- 
count the dispersion of the data in each radial bin. In Figure [4] 
we display for clarity only the dispersions for the total gas sur- 
face density distribution, which is of relevance for this paper. 
The CO-BIMA scale le ngth and its uncertainty are taken from 
Engargiol a et alj (|2003). The Hi surface density is remarkabl y 
constant with radius (as already noted bv lDevereux et al.ll 1 997b . 
with a value of ~10 M pc~ 2 , over most of the galaxy. Only be- 
yond 6 kpc, a cutoff is observed in the Westerbork data, likely 
due to interferometer losses where the ISM is more uniformly 
diffuse. The situation is different for the molecular gas: the CO 
emission decreases steeply outwards, implying a less efficient 
conversio n from atom i c to m olecular gas at larger radii. As dis- 
cussed bv lHever et al.l d2004) the molecular fraction is regulated 
by the balance between formation rate, related to the hydrostatic 
pressure in the galactic plane, and dissociation rate related to the 
ISRF. The hydrostatic pressure decreases rapidly outwards be- 
cause of the decrease of the stellar disk surface densit y, which 
domi nates over the gas surface density inside 5 kpc dCorbellil 
120031) . 

4. Dust in the disk 

4.1. UV and Optical extinctions 

To correct the UV and H a lumino s ities f or internal extinction we 
adopt the prescription bv lCalzettil d200ll) which relates Ay to the 
TIR and FUV luminosities : 



Ay = Cx 1.76 x log 



1 L(TIR) 
L68 X L(FUV) + 1 



(2) 



This gives an estimate of the average extinction of radiation from 
sources within the M 33 disk. The value of C depends on whether 
we consider extinction of the stellar continuum or of t he radia- 
tion fr om the ionised gas. In the model considered by ICalzettil 
( 200 1 ), the dust is driven to the outer part of star-forming regions 
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by first-generation massive stellar winds and supernova events. 
The dust will then produce a screen-like extinction for the radia- 
tion from central stars and from the ionised gas. Previous genera- 
tions of stars, still contributing substantially to the UV emission, 
will instead be interspersed within the dust, since they tend to 
diffuse out from the central region. Consequently, the extinction 
of the stellar continuum from older stars will be lower than the 
extinction suffered by younger stars and by the ionised gas. In 
practice, C — 1 for the ionised gas and C = 0.44 for the stel- 
lar continuum. In the case of M 33 this scenario is supported by 
the distinct morphology of the GALEX and Ha maps: the Hn 
regions are very localised, the UV much less so. However since 
about half of the Ha radiation in M 33 is diffuse we shall con- 
sider an average value: C = 0.7 for the ionised gas. This model 
assumes that the bulk of dust heating, which is re-radiated in the 
IR, comes from massive stars, the main contributors of UV pho- 
tons. As we will see later in this paper, this assumption might 
not be valid for the MIR radiation at 24 fim whose contribution 
to the TIR, however, is small, about 10%. 

Using a standard extincti on curve, Ay = 3.1 E(B - V), and 
following fSeibertetal.1 d2005l) . i.e. A FUV = 8.29 E(B - V), we 
can infer Apuv- Globally, Ay is found to be 0.30 mag (before ap- 
plying face-on correctio ns) and hence Apuy = 0.53 mag, consis- 
tent with the findings o f lBuat&Xul (fl996) which give an upper 
limit of 0.7 mag; hence the extinction corrected FUV luminosity 
of M33 is 6.3 x 10 42 erg s _1 . Using A NUV = 8.18 x E(B - V) 
dSeibert et al.l l2005). A NU y = 0.52 mag on average and the cor- 
rected NUV luminosity is 5.5 x 10 42 erg s _1 . 

In Fig. [5] we show the average face-on Ay value for the stel- 
lar continuum as a function of radius. Ay displays a rather con- 
stant value of about 0.27 mag up to 3 kpc, than drops gently to 
below 0.22 mag beyond 6 kpc. The drop of the UV luminosity at 
larger radii, due to the effective radial cut-off of the star forming 
disk, does not allow tracing the extinction further out, where in 
fact it may be considered negligible. 

Extinction is found to be globally rather low, in agree- 



begji 

ment with iTabatabaei et al.l d2007al) . Nevertheless, extinction 
can be enhanced locally: i n the cores of Hn regions values of 
Ay ~ 1 mag are common (Mass ev" et al J 1 19951 iDevereux et alj 
|1997|) . The optical extinction, measured from the Ha to WB ra- 
tio in selected H n regions, is found to have a large scatter with 
an average of 0.65 mag and with no clear radial dependence 
(Magrini et al. 2007b). Since, as noted above, the gas in H n re- 
gions should in principle suffer from a higher extinction than 
field stars, our average value of 0.25 mag for the stellar con- 
tinuum translates into Ay ^ 0.57 mag for th e H n regions, very 
close to the finding of Magrini etail d2007bl) . 



From the radial profile of Ay, we can recover the UV 
and Ha radial profiles corrected for extinction. These are 
shown in Fig. [3] and the relative scale lengths are given in 
Table [3] (FUV = 2.2 kpc ±0.1; NUV = 2.1 kpc ± 0.1; 
Ha = 1.8 kpc ± 0.1). As it can be seen, extinction corrections 
are small and do not dramatically affect the profile slope. Upon 
correction for extinction at each radial bin, we find the total FUV 
luminosity to be 7.0 x 10 42 erg s . We deem this value more ac- 
curate than the 6.3 x 10 42 erg s _I obtained above using the value 
of Apuv averaged over the whole disk. Similarly, the total NUV 
luminosit y radially correc ted for extinction is 6.3 x 10 42 erg s" 1 . 
Following Calzetti (2001) we retain a global estimate of extinc- 
tion for the ionized gas of 0.33 mag at Ha wavelength (C = 0.7). 
The corrected Ha luminosity is 4.2 x 10 ergs -1 . 
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Fig. 5. Ay face-on values for the stellar continuum averaged in 
radial bins of 0.24 kpc, as a function of galactocentric radius. 



4.2. Dust-to-gas ratio 

Since Ay is proportional to the dust column density, its ratio to 
the gas column density can be used to infer how the average dust- 
to-gas ratio varies across the M 33 disk. As discussed above, to- 
wards Hn regions extinction is higher than what Fig. [5] shows. 
Thus we expect the dust column density to be larger than what 
can be inferred from Figure|5] There is also the possibility of an 
additional dust component dispersed in the ISM and not heated 
by the young stellar population which has not been taken into 
account in the computation of Ay. 

In order to derive the dust column density from emission pro- 
cesses we shall derive the dust optical depth using the 160//m 
brightness and the dust temperature (see Section |4~3l for details 
on the dust temperature). The method used and the derived radial 
distribution of Ti6o, the optical depth at 160jum are very sim- 
ilar to those reported by Tabatab aei et alj d2007bl) for the same 
galaxy. The opacity decreases by a factor 10 from the center to 
the edge of the star-forming disk. We then convert the dust opac- 
ity to the global visual extinction through half of the face-on disk 
of M 33 using: 



1.086 

Ay = — - — x 2500 x T 16 o 



(3) 



Fig. [6] shows the ratio of Ay to the total gas surface density 
as a function of radius using the Ay values derived both from 
Eq. [2] and from Eq. [3] The gas is the sum of the H i (from the 
Arecibo data, except the first point which is from Westerbork) 
and molecular gas (from the FCRAO datafl There is a clear 
discrepancy: going radially outwards the dust-to-gas ratio de- 
rived from UV absorption does not show any clear dependence 
on radius. The dust-to-gas ratio derived from the IR emission de- 
creases steadily. Clearly both methods have a caveat: the use of 
Eq.|2]is valid under the assumption that all dust is heated by UV 
radiation, in the converse case it overestimates Ay. Also it de- 
pends on geometrical factors such as the clustering of hot stars, 
which emit the UV radiation. For the computation of Tigo we 



3 If only interferometric data were used to estimate the Hi density, 
due to the lower sensitivity and incomplete low spatial-frequency sam- 
pling, the dust-to-gas ratio would instead rise sharply in the outer part 
of the galaxy. 
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Fig. 6. The visual extinction normalized to the total gas column 
density as a function of the galactocentric radius. Plus symbols 
show this ratio for Ay derived from Eq. [2j open circles for Ay as 
from Eq.[3] 
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Fig. 7. Elliptically averaged 70 and 160 pun surface brightness 
and their ratio as a function of the galactocentric radius, follow- 
ing logarithmical bins. The two brightness profiles are in arbi- 
trary units and vertically displaced for clarity. 



have assumed a single dust component in each radial bin (ac- 
cording to the temperatures shown in Fig. [9j>: in the presence of 
several temperature components the characteristic temperature 
one derives is biased to the higher values, resulting in an under- 
estimate of the dust optical depth. Moreover the conversion fac- 
tor between the optical depth in the IR and the one in the visual is 
highly uncertain: values range from 2000 to 5000 for our Galaxy 
dWhitcomb et alJH98lHDraine & Led [l984). The dust composi- 
tion and size distribution in M 33 can also be different from that 
in our Galaxy and make precise calibration difficult (small par- 
ticules and highly reflective mater ials will tend to increase th e 
value of the conversion factor, e.g. lMivake & Nakagawalfl993h . 
Given these uncertainties, any conclusion about the average ra- 
dial variations of the dust-to-gas ratio and of the dust optical 
depth in M33 would be highly speculative. 
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Fig. 8. Elliptically averaged FIR emission per unit gas versus 
FIR colour. They are displayed for the 10 pm emission (upper 
panel) and for the 160 /im emission (bottom panel). In both pan- 
els, the galactocentric radius increases from the upper to the 
lower points. 



4.3. Dust temperature 

If the dust column density has a shallower radial decrease than 
the light profiles of the stellar population which is heating the 
dust, it is conceivable that the radial decrease of the 70 and 
160 pm emission is due to a decrease of the dust heating radi- 
ation field rather than to a decrease of the dust column density. 
If this picture is correct and the distribution of grain sizes does 
not change radially, any variation in the ratio of gas to 70 and 
160 pm emission is likely to be linked to dust temperature vari- 
ations and cooler grains are expected in the outer part of the 
galaxy. The radial behaviour of the 70/160 colour is shown in 
Fig.0 

There are clear oscillations, in coincidence with the spiral 
pattern, where the colour gets bluer due to the presence of large 
star-forming complexes. Overall this ratio decreases in the inner 
disk, where hot sources contribute substantially to the FIR emis- 
sion, then remains constant up to the edge of the star-forming 
disk. The ratio flattens out in the outer disk: in fact the scale 
length of the 160 pm emission, which was larger than that of 
the 70 pm emission in the inner disk, decreases and emission at 
both wavelengths falls off radially in the same way. At 7 kpc the 
data suggest an uprise but the large error bars and the possible 
disk distortion due to the warp prohibit any definitive conclu- 
sio n. By fitting an optic ally thin thermal spectrum (/q oc A~ 2 , 
see|Hirashita e t al.ll2007l and references therein) to the elliptical 
averages at 70 and 160 yum, we derive the integrated FIR emis- 
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Fig. 9. Dust temperature as a function of radius, estimated by 
fitting a blackbody to the 70 and 160 pm elliptical emissions, in 
bins of 0.24 kpc. 



sion by integrating an optically thin thermal spectrum between 
40 and 1100 pm, (which is the wavelength range which defines 
the FIR). The dust temperature, displayed in Fig. [9] decreases 
from 25 K in the innermost region to 21 K at about 4 kpc, then 
shows a flat plateau up to the edge of the star-forming disk where 
it increases. The integrated emission will be used in Sect. [6] to 
estimate the FIR luminosity and to derive the conversion factor 
between FIR luminosity and SFR. 

Figure [8] shows the ratio of the 70 or 160 pm surface bright- 
ness to the total gas column density as a function of the 70 to 
160 pm ratio. High values of the 70 or 160 pm surface bright- 
ness to gas column density ratio correspond to the inner disk. 
Here the strength of the emission at 70 and 160 pm is corre- 
lated primarily with dust heating and not with dust abundance 
since the FIR colour becomes redder as the dust emission per 
gas unit decreases. In the outer disk this is not so evident: the 
dust temperature stays constant out to the edge of the SF disk 
in the presence of a radially decreasing radiation field. If the 
dust column density does not change much and the UV radia- 
tion is the only heating source for grains, a change in the spec- 
tral energy distribution is required. But the lack of changes in 
the UV colour does not support this picture. A change in the 
dust abundance, as the steep radial decline of the 160/im emis- 
sion seems to suggest, is more likely, assuming that the dust size 
distribution does not vary. There is however some anomaly in 
the distribution of evolved stars which suggests that the build- 
ing up history of the M 33 disk has been rather complex and 
provides an alternative solution. An excess of Carbon Stars at 
large galactocentric radii has been foun d (e.g. lRowe et al.l i2005: 
iMcOuinn et al.ll2007t ICioni et al.l l2008). Accretion events in the 
past brought some metal poor material beyond 4 kpc thus favor- 
ing the production of Carbon Stars. The presence of additional 
radiation from those Carbon-rich giants in the outer disk can heat 
the dust in their proximity thus providing an additional contribu- 
tion to the 70pm surface brightness. We shall address better in 
the next Section the role of evolved stars in powering the unre- 
solved emission at shorter wavelengths. Some caution is needed 
here in drawing any definitive conclusion because at large radii, 
the 70 pm image has a high noise level due to residual artifacts 
(see Paper I). Concerning possible changes in the grain size dis- 
tribution we cannot fully constrain it with available data since we 
are not sensitive to emission from a very cold dust component. 



5. Discrete sources and diffuse emission 

In this Section, we derive the emission associated with discrete 
sources and compare this with the total emission in order to de- 
rive diffuse fractions as a function of radius. In general, M 33 is 
a galaxy with high diffuse fractions, the origin of which is still 
a subj ect of debate. The U V diffuse fraction of M33 was stud- 
ied bv lThilker et al.l d2005l) to obtain indications on the location 
and nature of the non-ionising photons. The FUV diffuse frac- 
tion appears to be remarkably constant (equal to about 0.65) as 
a function of radius. Only beyond the edge of the star-forming 
disk a slight increase of the FUV diffuse fraction could indi- 
cate that some of the emission was not produced locally but 
rather scattered by dust. The Ha diffuse fraction was derived by 
iHoopes & Walterbosl d2000l) and is about 0.4 throughout the SF 
disk. 

In | Paper II we used the SExtractor software 
dBerti n & Arn outsl 1 19961) to compile a catalogue of 515 
sources emitting at 24 pm. Following the process detailed in 
Paperjl we can extract sources at 8 pm, and in Ha. The cata- 
logue at 8 pm comprises 516 sources, and the Ha catalogue 413. 
The coarse resolution at 70 and 160 pm allows the extraction 
of only the most prominent sources, about 150 in total; in any 
case, the emission associated with discrete sources would be 
highly underestimated at these longer wavelengths and we will 
not discuss them further. 

Figure [10] illustrates the radial distributions of the integrated 
flux in discrete sources, in bins of 0.48 kpc. The upper panel 
shows that the IR flux at 8 and 24 pm contributed by dis- 
crete sources declines radially more steeply than the total (scale 
lengths of 0.86 + 0.06 and 0.88 ± 0.13 kpc for 8 and 24pm 
source number counts, respectively). This is not the case for the 
Ha sources, whose flux declines radially with a scale length of 
1.5 + 0.3 kpc, practically the same as the total Ha emission. 
Similar results are found for the contribution of UV sources. 
Therefore, it is only the IR flux associated with discrete sources 
which shows a steep radial decline. The very high IR diffuse 
fractions in the outer part of the SF disk are responsible for 
the smooth decline of the light at IR wavelengths; no edge is 
seen, compared to the Ha and UV distributions which drop more 
sharply beyond 6 kpc. 

The lower panel of Fig. [10] displays the results for the frac- 
tion of emission contributed by discrete sources with respect to 
the total. The Ha fraction shows some scatter but is nearly con- 
stant throughout th e disk and equal to 0.5 on average; this is 
in agreement with Thil ker et al.l d2005l ) who, however, found a 
slight increase towards the center. The fraction of IR emission 
in discrete sources is nearly the same for the two wavelengths, 
but in the inner disc (up to 3.5 kpc) nearly 40% of the emission 
is contributed by discrete sources, while at larger radii they can 
account for less than 20% 0. The drop of the IR flux associated 
with discrete sources is beyond 3 kpc, where spiral arms appear 
to lose their strength. In the absence of strong perturbations, the 
large dusty SF complexes disappear and smaller sources reside 
in less opaque clumps (see also Corbelli et al. 2008, in prepara- 
tion). 

5. 1 . Sources and diffuse TIR 

We will now derive the radial trend of the TIR over the entire 
disk and of a sample of selected sources on the 8 and 24pm 



4 The discrepant point at 24 jim just beyond 5 kpc is due to the con- 
tribution of IC 133. 
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Fig. 10. Upper panel: Infrared (in units of mJy pc~ 2 ) and Ha (in 
units of 10 13 erg cirT 2 s _I pc~ 2 ) radial surface brightness dis- 
tributions for discrete sources at 8 pm (green diamonds), 24 pm 
(yellow triangles) and Ha (black pluses). Bottom panel: Fraction 
of the total fluxes due to emission from discrete sources at 8 pm 
(green diamonds), at 24 pm (yellow triangles) and in the Ha line 
(black pluses) as a function of radius. 
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Fig. 11. The 24yum to TIR luminosity ratio as a function of the 
8 yum/24 pm flux ratio for 48 regions selected at 160 pm. The best 
linear fit to the data is shown by the plain (red) line. 



term in the above equation is small. This implies that the 8 pm 
emission from sources depends on their TIR emission but does 
not contribute appreciably to it, as expected when dust is not 
very hot and PAH features dominate the 8 pm band. Anyway, 
this result only holds for bright IR sources, those well visible, 
resolved and isolated in the 160/im map; the same might not 
apply to fainter sources. 

Using the radial profiles we can also derive the fit to the same 
relation but for the total emission (sources + diffuse component); 
we find : 



emissions. It is known that the relationship between 8, 24pm, 
and TIR luminosities depends on the proper ties of the galaxy 
considered and on the scale exam ined (e.g. lEngelbracht et alj 
120051: iPerez-Gonzalez et alj|2006l) . We will first examine if the 
TIR emission in M 33, associated with selected sources or with 
the overall disk emission, is correlated with the 8 pm emission. 
For this purpose we extracted sources from the 160yimi image 
and selected only those well fitted by closely circular PSF (ellip- 
ticity < 0.3). Then we measured them, within the aperture used 
at 160jum, also at 8, 24, and 10 pm. We finally used the formula 
in Eq.[T]to estimate the TIR emission of 48 sources. The 24pm 
emission from discrete sources amounts to about 8% of their to- 
tal TIR. Figure QT| shows the 24 pm to TIR luminosity ratio as a 
function of the 8 to 24 pm flux ratio. The best fit linear relation 
is : 

log F(TIR) = log F(24) + 1 .08 + 0.5 1 log ■ (4) 

The 0.51 + 0.06 coefficient is slightly lower than what 
ICalzetti et al.l (120051) derive in M 5 1 using a fixed and larger 
aperture size but similar to what lThilker et al.l (12007b) derive for 
M 33 using a sample of IR selected clumps on scales of 160 pc. 
Figure QT| shows that a correlation exists but also that the third 



log/^TIR) = log/ A (24) + 1.13 + 0.06 log(-^M . (5) 

\Iva{24)I 

The correlation between the TIR and the 8 pm emissions is 
in this case nearly absent, opposite to the trend found for se- 
lected IR sources. The absence of such correlation in the diffuse 
emission, that is away from bright H n complexes, is the result 
of a different radial trend of the 8 pm and of the IR emission at 
longer wavelength. 

Figure[T2] shows the remarkably good correlation, with slope 
about unity, between the TIR and the 24 pm fluxes when using 
elliptical averages. This does not necessarily imply that the pow- 
ering source in the two bands is the same. We shall see in the 
rest of this Section that if two stellar populations of different 
ages are responsible for heating the dust in the TIR and 24 pm 
band respectively, these results could imply that there has been 
no variation in the radial dependence of the SFR between the 
two epochs. 

5.2. Origin of the unresolved 24pm emission 

Discrete Hn regions in the innermost half of the SF disk are 
clearly powering at least 40% of dust radiation at 24 and 8 pm. 
Then, what is left must be associated with dust spread out in 
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Fig. 12. TIR versus the 24 /mi emission, from radial averages. 
The slope of the best fitting linear relation between their loga- 
rithmic values (continuous line) is 0.91. 

the disk and heated by the diffuse ISRF or it must be associ- 
ated with a "sea" of unresolved sources. The fraction of dust 
associated with H n regions decreases drastically with radius not 
only because the SFR per unit area is declining radially but also 
because they get dimmer in the IR beyond 3 kpc, as the drop 
in the bottom panel of Figure 9 shows. This raises the problem 
of the energy requirement at these wavelengths. Grains associ- 
ated with emission at 8 and 24 //m wavelengths are thought to 
be VSGs, not necessarily spatially coincident with the large-size 
grains responsible for 160 /mi emission. Additionally, at 8/im 
the emission might have a large contribution from PAHs. We 
therefore test two different hypotheses for powering the diffuse 
emission at 24 /urn: (a) the interstellar radiation which is giving 
most of the power by stochastic heating of VSGs; (b) evolved 
AGB stars, mostly carbon-type stars which are unresolved and 
have circumstellar dusty envelopes. 

The first issue is the location of any unresolved 24 /mi emis- 
sion. Radiation at this wavelength, in fact, must be produced 
in situ and cannot be imputed to scattering of radiation emit- 
ted elsewhere. This is easily verified by considering the ex- 
tremely low scattering cross section in t he MIR. For a stan- 
dard MR N dust mod el (Math is et alfl97 7p1 and the optical con- 
stants by lLaor & Draind d 1993b for the silicate material and by 
ZubkoetalJ $996) for amorphous carbon, we can estimate the 
scattering coefficient at 24 /mi to be about 1.4 x 1CT 6 times the 
extinction coefficient in the V band. As we have seen, the face- 
on average Ay value observed in the disk is about 0.25, which 
translates in a total optical thickne ss of the disk, after correcting 
for scattering (Bia nchi et al.l20 05). Ty — 0.35. This means that, a 
24 /mi photon may travel several Mpc in the M 33 disk before be- 
ing scattered into our line of sight by dust. We note, incidentally, 
that the values observed for Ay and for the total (H 1+H2) gas col- 
umn density Nh imply, for the dust model cited above, a mass 

s n(a) oc a~ 3 - 5 , a mi „ = 0.001 fim and a max = 0.3/mi, silicate and car- 
bonaceous grains in equal number. 



density ratio pdustlpH = 3.9 x 10~ 3 , slightly lower tha n in our 
Galax y where estimates span from 4 x 10~ 3 to 9 x 10 -3 dWhittetl 
1992 ), and in overall agreement with previous estimates, see 
llssa et al ] (119901) and references therein. 

5.2.1 . Emission of Very Small Grains 

In the diffuse ISM, away from circumstellar environments, 
stochastic, or single-photon, emission by VSGs is thought 
to be the main MIR continuum-emission mechanism (e.g. 
iDraine & Lil2007l) . Precise computations are not easy to perform 
and in what follows we illustrate a toy model aimed at a rough 
estimate of the VSG contribution to the diffuse 24 /mi emission 
inM33. 

The emissivity of a single grain when integrating over a time 
t long enough to smooth out the impulsive behaviour will be 

na 2 Q a v hs r T 

e v = — B v (T(t))dt (6) 

T Jo 

where Q" bs is the absorption efficiency of the grain and B v the 
Planck function at the grain temperature. We will model the 
temperature behaviour with time T(t) of a certain grain as a 
stochastic sequence of identical, non-overlapping, flat-top events 
of height T g and duration At, and T = otherwise. In each ab- 
sorption event the grain will attain a temperature T g which de- 
pends on the energy of the incoming photon and on the caloric 
properties of the grain. Given the photon, independently of the 
precise specific heat, the main factor will be the total number 
of atoms in the lattice since T g oc aT & , where 6 is the dimen- 
sionality. Small grains are not expected to be spheres, they may 
be flat or plain needles or, more probably, a mixture of shapes 
and dimensions: we assume an average dimensionality 6-2, 
i.e., planar. The grain will remain at temperature T g for the time 
At it takes to reradiate the absorbed energy. For a given incom- 
ing photon energy, na 2 (Q c f^) T g At = const., where (Q"^") is an 
average IR efficiency. For the same illuminating field, the num- 
ber of events per grain during t » At is N ev oc ma 2 (Qyy), 
where (g*) is an average efficiency for the absorption of the 

heating photons. For small grains of radius a < 100 A, we will 
always be in the Rayleigh regime and the absorption efficiency 
Qabs K a rp ne characteristic grain temperatures we are interested 
in are those for which AF x peaks at 24 //m, that is T g 150 K. 
Around this temperature, the Planck function at 24 jjm can be 
approximated, as a power law of T by B2a(T) oc T 3 9 . Putting 
all this together into Eq.[6] it is easily derived that the 24 /mi 
emissivity of the single grain varies as 624 x a 3 ' 2 - This is in re- 
markably good agreement with simulations of the VSG emission 
in the range 10 to 50 A of radius (Kruegel 2003). Obviously this 
relation will hold only in the range of grain radii which can be 
heated up to ~ 150 K by single-photon absorptions. The grains to 
be taken in consideration when heated by photons in the range 
900 < A < 5000 A are those containing roughly between 100 
and 600 atoms, which translates into 15 < a < 40 A. 

Following this line, upon averaging in time and composi- 
tion, we estimate the emissivity of a a = 40 A grain in the local 
ISRF to be 624 = 4.3 x 10~ 31 erg s _1 Hz _I sr _I . If we normalize 
the grain emission to this value and then integrate on the MRN 
size distribution between 10 and 50 A we can derive a volume 
emissivity and, given the observed gas-to-dust ratio, obtain the 
surface brightness at 24 /mi contributed by VSGs which, in this 
scheme, will be simply proportional to the gas column density 
and to the ISRF: cr™ c = 6.0 x lO~ 6 a H U mJy pc~ 2 , where cr H is 
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the H-gas column density in M pc~ 2 , and U the ratio between 
the actual ISRF and the one in the solar neighborhood 

In order to estimate the factor U, we derive in AppendixlAlan 
approximate solution for a disk-like atmosphere which allows us 
to derive the field at the midplane, given the surface brightness 
seen by an external observer. We will use this method to estimate 
the M 33 ISRF at vari ous wavelengths and compa re it to the local 
values in the Galaxy ( We ingartner & Dr aine 2001). 

Since the 3.6 /mi waveband is thought to directly measure the 
evolved stars photospheric emission, we may use it as a proxy 
of the overall stellar content and therefore of the ISRF itself. 
Given the stellar population responsible for the 3.6 //m emission, 
we will postulate a ratio between dust and stars scale heights 
P = 1/3; in M33 at this wavelength we estimate the optical 
half-thickness to/// - 0.02. With these assumptions, at 4kpc in 
the outer disk of M33, the observed 3.6 /mi surface brightness 
implies a midplane ISRF which is 0.31 times the local Galactic 
value. Alternatively, since the VSG stochastic emission is pow- 
ered mainly by UV photons, we may directly estimate the ISRF 
from the GALEX data. We assume f3 = 1/1 and to/// - 0.9 for 
both the NUV and FUV emission and in both cases we consis- 
tently derive, at 4 kpc, an ISRF 0.50 times the one in the solar 
neighborhood, in overall agreement with the 3.6 //m result. If we 
then assume that the ISRF in the outer M33 disk is about half 
the local Galactic value, we find that the predicted VSGs emis- 
sion falls short by a factor > 10, of the diffuse 24 /mi flux actu- 
ally measurecQ. Although exemplified here at a fiducial radius of 
4 kpc, this happens to be true at all radii, as well as for the total, 
unresolved 24 //m luminosity. We have verified that this result 
depends weakly on the limits used for the range of grain radii 
and on the slope of grain size distribution. Also, the VSG emis- 
sion of spherical grains will be lower than what we have here 
estimated for flat (6 = 2) geometries. 

5.2.2. Carbon stars 

To better understand the origin of the diffuse 24 /mi emission, 
we focus on a class of intermediate-age stars, the carbon stars. 
Carbon-rich giants, or carbon stars (CSs) are stars of low to in- 
termediate initial mass, 0.8 to 4 M with extended dusty circum- 
stellar envelopes mainly fo und in their asymptotic giant branch 
(AGB) evolutionary phase dHerwid2005l) . 460 carbon stars have 
been observed using the Spitzer infrared spec trograph (IRS) in 
the Magellanic Clouds bv lGroenewegen et al.l d2007l) and half of 
them in the LMC, an environment not dissimilar from the M33 
disk. While their AF ^ peaks around 8 /mi, they still show rather 
strong emission at longer wavelengths, with a AF^ at 100 /mi 
comparable to the one in the J bandpass. The observed spread 
in luminosity covers about 1 .7 dex with the most luminous ones 
being the reddest, as expected in the case of AGB stars; when set 
at the distance of M 33, the average 24 /mi flux of a single CS is 
F24 = 0.21 mJy. 

In order to estimate the possible contribution to the MIR 
emission by CSs we will mak e use of the Padova evolutionary 
tracks by Girardi et al. (2000) for a metallic content Z=0.00C 



6 In the same units iDraine & L.1 j2007h obtain cr^ G 2.0 x 
10^V«(/mJy per 2 . 

7 or a factor > 3 when adopting the more generous emissivities by 
IDraine & ll ll2007h . 

8 The tracks do not cover the full AGB phase but only up to the onset 
of thermal pulses. Under this respect, our estimates constitute a lower 
limit to the CS contribution; due to the short duration of the pulsing 
phase, the correction is of little importance. 



It is during the shell helium-burning phase that stars climb the 
red supergiant branch in the HR diagram by increasing their lu- 
minosity and lowering the effective temperature. The duration of 
this AGB phase varies with the initial mass, from about 10 7 to 
10 6 yr in the CS mass range. Similarly, the luminosity attained 
on the AGB depends on mass but three quantities are roughly 
constant with mass: (1) the effective temperature range from 
Tint — 4500 to T f ln si 3800 K; (2) the ratio between final and ini- 
tial luminosity Lfj n /Li n j 13 (that is, in log-log, different masses 
evolve along identical, vertically displaced, parallel tracks); and 
(3) the bolometric energy release J ACB Ldt =2 4.5 x 10 9 L yr. 

If L(m, £4) is the bolometric luminosity of a CS of initial 
mass m after a time t A since the entrance in AGB phase, we can 
define a time-averaged luminosity at mass m: 



{Urn)) 



J Q T ' ,<m) L(m,t A )dt A 
T A (m) 



(7) 



where T a (m) is the AGB phase duration for mass m. Given an 
Initial Mass Function, we can derive an overall mean luminosity 
for all CSs: 



Jo.8M G °< L ( OT » OT ~ gQfOT 
f 4.0M o 

,.„.. m "dm 

J 0.8 M 



(8) 



where a is the exponent of the IM F power law. H ereafter we 
adopt a = 2.3 for .5 < m ^ 100 M B dSalpeteril 19551) and 1.3 for 
0.1 < m < 0.5 M Q dKroupall200H) . For this IMF and the Padova 
tracks, we obtain ((L)) = 580 L . 

As we said, based on the observations in LMC, we know 
the average expected 24 /mi flux from a CS in M33, and we 
may think of this as a mass and time average over the whole 
CS population, as outlined above for the luminosity, that is 
((F24)) = 0.21 mJy. We also noted that the effective tempera- 
ture range of CS is rather restricted and quite similar for all the 
relevant masses. Given this circumstance, we may assume that 
with a good approximation ((F04)) oc ((L)) with a proportional- 
ity constant ({K 2 a-l)) = «^24»/«£» = 3.65 x lO^taTyLo" 1 . 

At this point, if F^itn, t A ) is the flux from the single CS, 
from a given region we expect a flux 



-4.0M G I r*T k (m) 

£24 = I N(m) < I FiA^m, t A ) dt A \ dm 
0.8 m [Jo 



(9) 



where N(m) is the differential star-formation rate in that re- 
gion, admitting it is constant with time (both integrally and 
in mass distribution). If we now assume that F^irn, t A ) - 
{{K24-L))L{m,t A ), the theoretical tracks can be used to evalu- 
ate Eq. [9] and we obtain £ 2 4 = 2.05 x 10 5 M7-mJy, where M T 
is the total (from 0.08 to 100 M ) SFR in the region of interest 
in M yr _I . Note that, despite the rather strong assumption of 
proportionality between F2\{m, t A ) and L (m, t A ), the final result 
relies on estimates of J Ldt along the AGB phase which, as said 
above, are quite robust and constant. 

It turns out that the SFR required to explain the 24 /mi un- 
resolved emission in M 33 with the sole contribution of CSs is 
quite compatible with what is computed by the usual SFR indi- 
cators, as will be seen in the following. Note also that the most 
luminous CSs could well have been counted among the faint dis- 
crete sources, and that a small fraction o f carbon stars might not 
have dusty shells (Mc Ouinn et al.ll2~007l) . As a result the 24 /mi 
unresolved emission due to CS should give a lower SFR than the 



12 



S. Verley et al.: Star formation in M 33: multiwavelength signatures across the disk 



effective one. The total unresolved 24 pm flux of M 33, ~ 35 Jy, 
translates into a global SFR of 0.17 M yr~'; the surface bright- 
ness in the outer disk at 4kpc is ~ 4 x 10~ 4 mJy pc~ 2 requiring 
M T ^2.0M o Gyr" 1 pc- 2 . 

6. Star formation rates 

How well tracers at different wavelengths measure the actual 
SFR has b een examined by a number of papers based on Spitzer 
data (e.g. ICaizetti et al J 120051: iKennicutt et alj|2007l for M 51). 
In this Section, we compare the SFRs derived from emission at 
different wavelengths, as a function of the galactocentric radius 
(in bins of 0.24 kpc) and integrated over the whole SF disk. To 
compute the SFR from Ha and UV luminosities corrected for 
extinction (see Sect. [3) we use the stellar p opulation synthesis 
model in Starburst 99 dLeifherer et al J [19991) . for a continuous 
star-formation model with solar metallicity and the IMF cited in 
the previous Section. 

6.1. SFR from Ha and UV wavelengths 

The conversion between Ha luminosities and SFR given by the 
stellar population model is 

SFR(Hff) [M® yr~'] = 8.3 x 10~ 42 L(Ha) [erg s" 1 ] . (10) 

This is within 5% of the value given bv iKennicutJ (1 19981) . 
showing that the SFR dependence on Ha is rather insensitive to 
the assumptions about the lower-mass form of the IMF. The ra- 
dial profile is shown in Fig.[13]and the total integrated SFR over 
the whole SF disk (R<1 kpc) is SFR(Ha) = 0.35 M yr _1 which 
correspond to the average SFR over about the last 10 Myr. To 
trace the SFR over a wider timescale, namely in the last 100 Myr, 
we convert the FUV luminosity according to the same popula- 
tion synthesis model at 1550 A as: 

SFR(FUV) [M yr _1 ] = 8.8 x 10~ 44 L(FUV) [erg s" 1 ] . (11) 

The above conversion (equivalent to 1.7 x 
10~ 28 L V ( FUV) [erg s~' H z" 1 !) is also similar to the one 
given bv lKennicuttl (Il998l) (although we are using a slightly dif- 
ferent IMF and wavelength range). The radial profile of the SFR 
derived from the FUV emission is shown in Fig.[T3]and the SFR 
integrated over the SF disk is SFR(FUV) = 0.55 M yr _1 , higher 
that the one derived from Ha. The SFR declines radially with a 
scale length of 2.2 + 0. 1 kpc in the FUV and 1 .8 ± 0. 1 kpc in Ha, 
a marginal difference considering the uncertainties. However 
there are two reasons why the two might indeed be different: 
one is because we are sampling two different SF epochs; the 
second is that Ha depends on the interstellar medium volume 
density and a radial decrease of this might suppress Ha an d 
also induce a sharper Ha SF edge (Elmegreen & Hunter 2006). 

6.2. SFR from IR wavelengths 

While in dusty environments IR-based SF tracers are quite reli- 
able, in M 33 this assumption is not obvious and requires cau- 
tion in interpretation. Although within the single SF region most 
of the dust emission is powered by young stars, the situation is 
more complex when we integrate over the whole disk which, as 
said before, has an important diffuse fraction. In M33 the emis- 
sion at 70 and 160yum can be powered by the UV in the diffuse 
ISRF or by single H n regions: both cases are related to the pres- 
ence of stars of young to moderate age. We expect, therefore, 



the FIR luminosity to be a valid SF tracer and yield a SFR sim- 
ilar to the one based on FUV data, once corrected for extinc- 
tion. If we compute the FIR luminosity by fitting an optically 
thin (k,i oc /T 2 ) thermal spectrum to the 70 and 160 /im elliptical 
averages and then integrate this between 40 and 11 00 yum (see 
Sect. 3J, we find that it correlates well with the TIR luminos- 
ity computed via Eq. Q] The average ratio between the TIR and 
the FIR luminosities is 1.40 + 0.04 across the whole disk. The 
FIR emission thus modeled, integrated over the disk (< 7 kpc), 
amounts to 4.6 x 10 42 erg s _1 . 

If we choose the extinction corrected FUV as a reference 
for the SFR, we can derive the conversion factors for the other 
wavebands. We did this by using a^ 2 minimization over the set 
of radial bins within 6 kpc, to avoid the most noisy data. For the 
FIR we obtain: 

SFR(FIR) [M e yr" 1 ] = 13 X 10~ 44 L(FIR) [erg s" 1 ] , (12) 

which is in the upper range of the values found by iBuat & Xul 
( 1996) for late-type galaxies based on IRAS and UV flux mea- 
surements. A relatively high conversion factor is indicative of a 
relatively low dust content. 

Using the same technique but for the TIR fluxes, we find that 
in order to match the FUV SFR: 

SFR(TIR) [M e yr" 1 ] = 8.8 x 10 44 L(TIR) [erg s" 1 ] . (13) 

As seen in Fig. [13] the FIR- and TIR-based SFR radial pro- 
files show a smoother drop at the SF disk edge possibly due 
to their origin from more diffuse components, dust, and general 
ISRF. 

If one is convinced that also the 24 pm is powered by young 
stars, then normalizing the SFR based on the 24 pm emission to 
the FUV SFR, one obtains: 

SFR(24) [M @ yr" 1 ] = 1.2 x 10~ 42 L(24) [erg s" 1 ] . (14) 

This factor is about three ti mes high er with respect to the re- 
lation given bv lWu et al J (120051) (see also lYun et al.l200ll) . which 
was derived by comparing IR and radio emission in the Spitzer 
First Look Survey Galaxies. They find, in fact: 

SFR(24) [M @ yr" 1 ] = 3.9 x 10~ 43 L(24) [erg s" 1 ] . (15) 

Using this latter relation the global SFR from 24 pm emis- 
sion in M 33 is only 0.2 M yr" 1 and the relative SFR per unit 
area is shown in Fig.[T3] Eq.[T5l clearly underestimate the SFR. 
In the inner regions the 24 pm and Ha SFRs follow the pattern of 
the galaxy: they are enhanced along the spiral arms and decrease 
monotonically away from the center. Beyond 4 kpc the 24 pm, 
emission fraction linked to discrete sources drops. Diffuse emis- 
sion dominates which, as shown in the previous Section, can 
be powered by more evolved stars and linked to the SFR about 
1 Gyr ago rather than to recent episodes of SF. This explains the 
absence of any drop at the SF disk edge due to stellar diffusion 
of the more evolved population. Hence, the 24 pm emission in 
late-type galaxies like M 33, with a few very large SF complexes 
and a large diffuse component, depends upon the presence of 
dusty envelopes around evolved stars. If the contribution in hot 
discrete sources is removed, it can give a lower limit to the SFR 
about 1 Gyr ago: the exact rate depends upon the fraction of 
evolved stars which have developed a dusty envelope. 
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Fig. 13. SFR per square parsec as a function of radius using 
Eq. QT| for converting into SFR the FUV luminosity (open cir- 
cles), Eq. Q3] for the 24 pm luminosity (open triangles) Eq. [10] 
for the Ha luminosity (filled diamonds) (upper panel), Eq. 1 1 
for the FIR luminosity (open squares) and Eq. [13] for the TIR 
luminosity (filled circles) (bottom panel). 



7. Summary and conclusions 

This second paper of the series focuses on various SFR diagnos- 
tics across the disk of the local spiral galaxy M33, and on the 
origin of the large IR diffuse emission. We investigate the radial 
trends of the IR emission at various wavelengths using Spitzer 
data, complementing these with GALEX data for the UV, with 
Ha data, and with data on the atomic and molecular gas distri- 
butions. Hereafter, we summarise our main results. 

- The Spitzer radial profiles at 3.6, 4.5, 24, and 70 pm are con- 
sistent with a single exponential scale length throughout the 
star-forming disk. This is ~ 1 .55 kpc for the (old) stellar pop- 
ulation emission and ~1.75 kpc for the dust emission at 24 
and 70 pm. Two different slopes for the inner and outer star- 
forming disk are found instead at 5.8, 8.0 and 160 pm. At 
these wavelengths the IR brightness declines more steeply 
beyond 3.5 kpc and drops at the edge of the star-forming 
disk. The scale lengths more directly linked to SF diagnos- 
tics (Ha, FUV and NUV) are about 2 kpc, comparable to the 
IR ones in the inner disk, but longer than these in the outer 
star-forming disk. This implies that the dust column density 
is decreasing beyond 3.5 kpc, a trend confirmed by our esti- 
mate of the extinction from the TIR to FUV ratio. The 8 pm 
emission is dominated by PAH features which show instead 
a faster decline with radius. 



Within the star-forming disk, that is within 7 kpc of radius, 
we estimate the TIR, FUV, NUV, and Ha luminosities to be 
6.8 x 10 42 , 7.0 x 10 42 , 6.3 x 10 42 , and 4.2 x 10 40 erg s" 1 re- 
spectively after correcting for extinction. Extinction correc- 
tions are small and reflect the rather low dust content of the 
galaxy (Ay ~ 0.25, for stellar continuum averaged in rings of 
width 0.24 kpc). Extinction can be enhanced locally (cores 
of Hn regions, for instance); however, the absence of large 
dusty complexes beyond 3 kpc, likely due to the weakness 
of the spiral pattern, implies that at large radii SF proceeds 
in smaller clumps with low extinction. Dust opacity derived 
from the 160/mi emission under the assumption of single 
temperature component has a steeper radial decline than dust 
opacity inferred via the extinction model. 
A substantial amount of the emission of the M 33 galaxy is in 
diffuse form at several wavelengths. The FUV diffuse frac- 
tion is equal to about 0.65 and does not depend on radius. 
The Ha diffuse fraction is more scattered but nearly constant 
throughout the disk of the galaxy and equals to 0.5 on aver- 
age. In the IR, at 8 and 24 pm, about 40% of the total emis- 
sion of the inner disk comes from discrete sources while this 
contribution drops to 20% in the outer disk. The dust tem- 
perature, estimated from the 70 and 160 pm emission, de- 
creases from 25 to 21 K from the center to 4 kpc and, after a 
plateau, it might increase at the edge of the star-forming disk. 
The average ratio between the TIR and FIR luminosities is 
1 .40 + 0.04 across the whole disk. 

We examine the SFR of M 33 from various tracers (extinc- 
tion corrected Ha and FUV, TIR, FIR, 24 pm). We estimate 
the total integrated SFR to be 0.45 + 0.1 M yr -1 . From the 
Ha luminosity (which accounts for SF in the last < 10 Myr) 
we infer a lower value, SFR(Ha) = 0.35 M yr -1 , while the 
FUV emission (which accounts for SF in the last < 100 Myr) 
gives SFR(FUV) = 0.55 M yr -1 ; the difference is small 
with respect to the uncertainties, so that we cannot conclude 
definitively that SFR is declining in the last 100 Myrs. If we 
calibrate the SFR from FIR and TIR emission through the 
values from UV data, we obtain conversion factors in the up- 
per range of what is generally found for late-type galaxies, 
likely due to the low dust content of M 33. 
A SFR diagnostic based on the 24 pm emission in M 33 gives 
a conversion factor more than three times higher than the 
standard ones. To estimate the validity of the 24pm emis- 
sion as a SFR tracer away from H n complexes, we inves- 
tigated the nature of its high diffuse fraction. We find that 
VSGs emission, as powered by the diffuse ISRF, accounts for 
only 10% of the diffuse 24 pm emission, while unresolved, 
AGB stars of intermediate mass, Carbon Stars, through their 
dusty circumstellar envelopes, can provide most of the dif- 
fuse 24 pm emission. If indeed evolved stars contribute sub- 
stantially to the 24/i/m brightness, emission at 24 pm can 
trace SFR about 1 Gyr ago. In this paper we only give a 
lower limit to it of about 0.20 M yr . This estimate takes 
into account only the diffuse 24 pm emission and does not 
consider the most luminous CSs, which might be resolved 
as discrete sources, as well as the fraction of CSs without 
dusty shells. With the help of future surveys of evolved stars 
in M 33, it should soon be possible to complement our study 
and obtain a better view of the time evolution of the SFR in 
the Triangulum galaxy. 
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Table A.l. Midplane to emergent intensity ratio J V (Q) /I(0,fi). 
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Appendix A: A simplified exponential atmosphere 

We model the disk as a plane parallel medium symmetric respect 
to midplane. In case of isotropic extinction and emissivity, at 
vertical distance z from midplane, the transfer of the intensity 
I v {z, 9) along a direction at angle 9 from the vertical will follow 



dl v (z,0) 
Xv (z) dz 



Sy(z)-I V (Z,ff) 



(A.l) 



as usual fi - cos 9, Xv is the volume extinction coefficient, and 
S v is the source function; S v (z) = e v (z)/Xv (z)> where e v (z) is the 
volume emissivity. If r v is the vertical optical thickness r v (z) 
Xv (z') dz' , then 



V -J = I V {T,H) -S V {T) 

dT 



(A.2) 



If we assume, though not strictly required, that the intensity at 
midplane is also isotropic, I v (z = 0,6) = J V (Q), then the intensity 
emerging from the atmosphere (z = +°o, t = 0, ft > 0) will be 



I Y (P,H) = J Y (0)e 



■To//' 



' T ° At 



(A.3) 



where to = t(z = 0) is the optical half-thickness of the disk. At 
the same time, at midplane (z = 0, r = to, ft < 0), if there is no 
incoming radiation from outside the disk, we have 



7 v (T ,//)=/y(0) 



X 



T ° At 

S(T)e- (T °- T)/ » — 



(A.4) 



(A.5) 



If the extinction and emissivity volume coefficients depend 
exponentially on z, so will be S v (z) and the integrals in 
Eq.s lA.3IIA.4l can be evaluated in terms of Incomplete Gamma 
functions y (/?, + To///), where f3 is the ratio between the extinc- 
tion and emissivity scale heights,/? = h ext /h em . By comparing 
the two solutions we find the emergent intensity 7 v (t = 0, fi) in 
terms of the one at midplane J V (Q), and vice versa: 

{-Yfy{P, -To/fx) 

Some numerical values of 7 V (0) //(0,/i), the ratio between mid- 
plane and emergent intensity are shown in Table lA. ll For /3 = 1 
we recover the values for the particular case of constant 5,,. 
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